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ULTRAFINE  MICRO STRUCTURE  COMPOSITES 
PREPARED  BY  CHEMICAL  VAPOR  DEPOSITION 


I.  V  SUMMARY 

^  It  is  our  goal  to  develop  an  analytical  model  which  predicts 
the  experimental  conditions  which  permit  the  preparation  of 
ultrafine  microstructure  ceramic  composites  by  chemical  vapor 
deposition  (CVD).  The  model  is  to  be  based  on  classical  thermo- 
dynamics,  mass  transport,  kinetic,  nucleation,  and  growth  theory. 
The  model  will  be  validated  for  two  dispersed  phase  composite 
systems.  One  system  will  be  boron  nitride  plus  aluminum  nitride. 
The  other  will  be  hafnium  diboride  plus  silicon  carbide.  The  BN 
+  AlN  system  is  of  interest  for  radomes,  windows,  and  tribolo¬ 
gical  applications  while  HfB2  +  Sic  offers  potential  as  a  high 
temperature  oxidation  resistant  material  including  oxidation 
protective  coatings  for  carbon-carbon  composites T~) 

A  thorough  literature  search  was  completed  for  chemical 
vapor  deposition  of  BN;  AlN,  and  codepoeltion  of  dispersed  phase 
composites  including  HfB2  +  SiC.  The  literature  on  modeling  of 
the  chemical  vapor  deposition  process  was  reviewed.  Existing 
theory  applicable  to  our  model  is  reviewed  in  this  annual  report. 

Thermodynamic  analyses  of  the  BN  +  AIN  and  the  HfB2  +  SiC 
systems  have  been  completed.  These  computerized  studies  indicate 
that  both  composite  systems  can  be  prepared  by  codeposition  using 
commercially  available,  conventional  reagents*^ Results  for  the 
HfB2  +  SiC  system  were  published.  '  \ 

Mass  flowmeters  and  a  pressure  control  system  were  pur¬ 
chased  and  installed  on  the  Georgia  Tech  CVD  System.  This 
equipment  permits  accurate  control  of  reagent  flow  rates  and 
total  gas  pressure.  Deposition  experiments  have  emphasized  BN, 
AIN,  and  codeposition  of  these  two  materials  to  form  a  dispersed 
phase  composite.  Preliminary  analysis  of  the  codeposited 
material  by  ESCA  reveals  the  desired  chemical  constituents  (B, 
Al,  ana  N)  but  more  detailed  characterization  (XRD,  SEM,  and 
electron  microprobe)  is  required  to  verify  the  existence  of  the 
desired  dispersed  phase  composite. 


II.  INTRODUCTION 

Because  the  properties  of  composites  are  frequently  superior 
to  those  obtainable  with  single  phase  materials,  it  is  of 
interest  to  explore  fabrication  of  composites  by  chemical  vapor 
deposition.  Multiphase  composites  have  been  prepared  by  chemical 
vapor  deposition  by  either  codeposition,  infiltration  of  fibrous 
or  particulate  preforms,  or  by  alternately  depositing  one  or  more 
materials.  Our  work  addresses  codeposition  of  ceramic  compos¬ 
ites,  i.e.,  simultaneous  deposition  of  multiple  phases.  In 
addition  to  the  matrix  phase,  one  or  more  dispersed  phases  are 
present.  The  shape  of  the  dispersed  phase(s)  may  be  spherical, 
platelike,  rodlike  (whiskers),  or  otherwise.  We  have  recently 
reviewed,  in  an  open  literature  publication [ 1 ) ,  codeposition  of 
dispersed  phase  composites  by  chemical  vapor  deposition. 

Several  reasons  exist  for  the  current  high  interest  in 
composites.  For  ceramics,  the  major  thrust  is  to  reduce  brittle¬ 
ness,  i.e.  increase  fracture  toughness.  Extensive  progress  has 
been  made  in  improving  mechanical  properties  of  particulate, 
whisker,  and  fiber  reinforced  ceramics  prepared  by  chemical  vapor 
deposition  and  by  techniques  other  than  chemical  vapor  deposi¬ 
tion,  e.g.,  powder  blending  followed  by  sintering  or  hot  press¬ 
ing.  In  addition,  silicon  carbide  whiskers  or  metal  particulates 
have  been  added  to  a  wide  variety  of  materials. 

Composites  which  contain  a  lubricating  phase  as  well  as  a 
hard  phase  (e.g.  BN  +  AIN)  offer  the  promise  of  tailoring 
properties  such  as  friction  coefficient,  hardness,  and  resistance 
to  wear  and  erosion.  Addition  of  a  hard  phase  to  zinc  sulphide 
or  zinc  selenide  would  also  be  expected  to  improve  the  wear  and 
erosion  resistance  of  these  important,  but  soft,  window  mate¬ 
rials.  In  addition  to  improved  tribological  properties,  enhanced 
resistance  to  oxidation  has  also  been  demonstrated  for  HfB2  + 
SiC.  This  dispersed  phase  composite  has  greater  resistance  to 
oxidation  than  either  HfB2  or  SiC  alone. [1]  Electrical,  mag¬ 
netic,  optical,  and  thermal  properties  can  also  be  favorably 
influenced  by  selection  of  the  proper  composite  system.  Deposi¬ 
tion  of  composite  coatings  .is  also  of  interest  because  of  the 
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possibility  of  tailoring  the  coatings  thermal  expansion  coeffi¬ 
cient.  Increased  coating  toughness  and  strength  and  control  of 
the  thermal  expansion  mismatch  between  the  substrate  and  coating 
should  contribute  to  better  coating  adherence. 

A  major  difficulty  in  conventional  preparation  of  dispersed 
phase  ceramic  or  metallic  based  composites  is  achieving  the 
desired  microstructure,  i.e.,  numerous  uniformly  dispersed  small 
dispersoid  particles.  This  type  of  structure  is  difficult  to 
obtain  by  melting  or  by  powder  blending  followed  by  consolidation 
since  the  dispersed  particles  tend  to  be  too  large  and  they 
segregate.  A  further  problem,  in  the  case  of  powder  mixtures,  is 
that  they  tend  to  be  difficult  to  sinter  to  acceptably  high 
densities  particularly  for  covalently  bonded  compounds,  such  as 
SiC  and  and  for  refractory  carbides  and  nitrides  and 
refractory  metals  in  general.  One  method  for  circumventing  these 
problems  is  to  utilize  heat  treatment  to  cause  precipitation  of 
the  desired  dispersed  phase.  While  this  is  practical  for  several 
metals,  unfortunately  not  many  ceramic  systems  exhibit  the 
required  solid  solubility.  Chemical  vapor  deposition  offers  a 
very  versatile  process  which  should  be  applicable  to  a  large 
number  of  matrix-dispersoid  chemical  combinations.  Thus, 
developing  a  basic  understanding  of  the  interrelationships 
between  processing,  structure,  and  properties  of  dispersed  phase 
composites  prepared  by  chemical  vapor  deposition  is  a  worthy 
objective. 

Obviously  if  the  goal  is  tc  deposit  a  multiphase  composite, 
the  process  conditions  must  be  selected  and  controlled  in  order 
to  permit  the  deposition  of  the  desired  phases.  As  described 
below,  thermodynamic  calculations  are  extremely  valuable  for  this 
purpose.  Reaction  kinetics,  mass  transport,  and  crystal  nuclea- 
tion  and  growth  are  also  very  important  since  they  influence  the 
microstructure  and  therefore  most  properties  of  the  deposit. 
Important  microstructural  features  for  the  matrix  and  dispersed 
phase  are  grain  size,  shape,  and  orientation.  These  attributes 
can  be  deliberately  altered  by  variation  of  the  process  condi¬ 
tions  most  responsible  for  influencing  reaction  kinetics  and 


crystal  nucleation  and  growth.  Progress  in  modeling  the  CVD  of 
BN  +  AIN  and  HfB2  +  Sic  is  described  in  the  remainder  of  this 
report  along  with  a  description  of  our  experimental  results  on 
deposition  of  BN,  AIN,  and  BN  +  AIN. 

III.  THERMODYNAMICS 

The  extensive  understanding  of  the  thermodynamics  of 
chemical  vapor  deposition [2-6]  is  of  great  value  in  predicting 
the  process  conditions  which  should  lead  to  codeposition  of  the 
desired  solid  phases.  Computerized  calculations  predict  the 
equilibrium  phases  based  on  minimization  of  the  free  energy  of 
all  possible  gaseous,  liquid,  and  solid  species. 

We  have  performed  numerous  thermodynamic  calculations  for 
the  BN  +  AlN  and  HfB2  +  SiC  systems.  The  calculations  were 
performed  usi..g  the  SOLGASMIX-PV  program[7]  which  we  modified  in 
order  to  permit  its  use  on  a  personal  computer.  Results  for  both 
chemical  systems  are  summarized  below. 

1.  BN.+  AIM 

The  goal  of  the  thermodynamic  calculations  is  to 
predict  reagent  concentrations,  pressures,  and  temperatures  which 
will  permit  codeposition  of  the  desired  BN  +  AIN  composite. 
Also,  the  information  learned  about  gaseous  species  assists 
understanding  and  modeling  of  chemical  kinetic  and  mass  trans¬ 
port.  Composites  covering  a  range  of  AlN  contents  are  desirable. 
For  example,  addition  of  a  small  amount  of  AIN  to  BN  should 
improve  its  erosion  resistance  while  additions  of  BN  to  AlN 
should  result  in  better  tribological  properties  because  of  the 
lubricating  nature  of  BN. 

The  thermodynamic  calculations  were  performed  for  the  BCI3- 
AICI3-NH3  reactant  system  over  the  temperature  range  1000  to 
1700  K.  The  results  were  very  encouraging  in  that  conditions 
were  found  for  which  it  is  predicted  that  BN  +  AIN  can  be 
codeposited.  The  calculations  show  that  it  should  be  possible  to 
prepare  BN  +  AIN  composites  which  vary  in  AlN  content  from  0  to 
100%. 


A  draft  report  describing  the  results  of  the  BN  +  AIN 
thermodynamic  calculations  has  been  prepared.  The  type  of 
information  learned  from  this  study  is  illustrated  by  the 
deposition  diagram  of  Figure  1.  This  diagram,  which  resembles, 
but  differs  from,  a  traditional  phase  diagram  shows  the  equilib- 
rium  solid  phase(s)  as  a  function  of  reagent  concentration.  The 
region  of  importance  is  the  two  phase  region  identified  as  BN  + 
AlN  in  Figure  1,  located  near  the  right  side  of  the  composition 
triangle.  Within  this  region  codeposition  of  BN  +  AIN  is 
predicted.  The  composition  of  the  composite  varies  from  0  to 
100%  AlN  on  moving  from  left  to  right  across  the  two  phase 
region. 

The  deposition  diagram  in  Figure  1  is  for  a  specific 
temperature  and  total  gas  pressure  (1400  K  and  1.0  atm)  for  the 
reagent  system  BCI3-AICI3-NH3 .  Numerous  additional  diagrams  were 
constructed  in  order  to  evaluate  the  influence  of  temperature, 
pressure,  and  hydrogen  concentration  on  the  codeposition  process. 
For  example,  the  influence  of  temperature  and  the  ratio  of  A1  to 
B  in  the  input  reactants  is  given  by  Figure  2.  This  graph  will 
be  useful  in  selecting  process  conditions  appropriate  for 
deposition  of  BN  +  AlN  composites  of  specific  BN  contents. 
Codeposition  appears  to  be  feasible  over  a  rather  wide  range  of 
temperatures  (1000-1400  K)  and  pressures  (0.1-1  atm).  The 
calculations  also  revealed  that  the  deposition  efficiency  was 
higher  for  BN  than  AIN  and  that  addition  of  hydrogen  increased 
the  deposition  efficiency  for  AlN.  Additionally,  the  influence 
of  partially  or  completely  replacing  BCI3  with  B2H5  was  analyzed. 
It  was  learned  that  both  BCI3  and  B2H5,  or  mixtures,  should  be 
appropriate  for  codeposition  of  BN  +  AlN. 

The  thermodynamic  calculations  are  also  valuable  in  that 
they  give  the  partial  pressures  of  the  numerous  gaseous  species 
present  during  the  deposition  process .  An  illustrative  plot  of 
pressure  versus  reactor  temperature  for  the  major  gaseous  species 
appears  in  Figure  3.  Such  plots  will  be  helpful  in  future 
kinetic  and  transport  modeling  studies  and  have  already  been 
useful  in  helping  us  to  understand  the  codeposition  process.  For 


example,  the  rather  high  stability  (partial  pressure)  of  the 
AICI3,  AlCl2and  A1C1  species  accounts  for  the  lower  deposition 
efficiency  of  AIN  compared  to  BN.  This  gave  us  the  clue  that 
addition  of  hydrogen  and/or  NH3  to  the  normal  reactant  stream 
would  lead  to  higher  AIN  deposition  efficiencies  since  their 
presence  would  increase  the  stability  of  HC1  and  decrease  that  of 
the  A1C1X  species.  This  happens  because  the  source  of  the 
chlorine  for  formation  of  HC1  is  the  A1C1X  species.  Subsequent 
thermodynamic  calculations  verified  this  hypothesis. 

2.  HfB2  +  SiC 

The  dispersed  phase  composite  HfB2  containing  26  mole 
percent  (20  vol.%)  SiC  has  been  shown  to  have  superior  resistance 
to  high  temperature  oxidation  and  thermal  shock  than  either  HfB2 
or  SiC  alone.  [1]  For  this  reason  and  also  because  the  borides 
have  superior  resistance  to  thermal  shock  compared  to  most 
ceramics,  the  composite  HfB2  +  SiC  is  of  interest  as  a  high 
temperature  structural  material. [8-11]  The  material  has  been 
prepared  by  powder  blending  followed  by  hot  pressing  but  attempts 
to  prepare  coatings  of  HfB2  +  SiC,  or  of  the  analogous  material 
ZrB2  +  SiC,  by  chemical  vapor  deposition  Mve  been  unsuccess¬ 
ful.  [10,12]  The  deposits  did  not  consist  of  the  desired  phases. 
Fenter[10]  reports  that  a  dense  ZrB2  matrix  containing  B4C, 
instead  of  SiC,  was  obtained.  The  reagents  used  in  that  work 
were  ZrCl4,  BCI3,  SiCl4,  and  H2.  The  deposition  temperature  was 
not  given. 

In  an  attempt  to  learn  how  to  simultaneously  deposit  the 
desired  HfB2  +  Sic  composite  by  CVD,  we  have  studied  the  thermo¬ 
dynamics  of  the  process.  Thermodynamic  calculations  were 
performed  for  the  reactant  system  HfCl4,  BCI3,  and  CH3SiCl3 
(methyltrichlorosilane  or  MTS)  again  using  the  SOLGASMIX-PV 
program.  A  quantity  of  hydrogen  was  added  equal  to  ten  times  the 
sum  of  the  number  of  moles  of  HfCl4,  BCI3,  and  MTS  present. 

The  results  of  the  calculations  are  very  encouraging  in 
that,  barring  kinetic  intervention,  several  conditions  were  found 
for  which  it  is  predicted  that  the  desired  HfB2  +  Sic  composite 
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will  be  deposited.  For  example,  the  CVD  deposition  diagram  in 
Pigura  4  shows  the  solid  phases  dapoaitad  for  tha  HfCl4-BCl3-MTS 
raactant  syatem  for  a  tamparatura  of  1200  K  and  a  total  eyatam 
praaaura  of  1  atm.  Tha  desired  ragion  whara  HfB2  and  Sic 
codapoait  ia  locatad  near  tha  MTS  cornar  of  tha  diagram.  Tha 
phase  boundaries  in  Figure  4  ware  determined  by  performing  tha 
thermodynamic  calculation  for  each  of  231  points  locatad  within 
tha  composition  triangle  representing  all  possible  5  mole  % 
increments  for  tha  reactants.  In  order  to  more  accurately  define 
tha  location  of  tha  phase  boundaries ,  calculations  ware  performed 
at  2  mole  %  increments  for  tha  MTS  rich  cornar  of  the  diagram. 
The  resulting  deposition  diagram  appears  in  Figure  5.  Tha 
contour  lines  within  the  HfB2  +  SiC  region  of  Figure  5  give  tha 
mole  %  SiC  predicted  for  the  two  phase  composite.  They  show  that 
it  should  be  possible  to  prepare  deposits  having  a  broad  range  of 
compositions,  including  the  desired  composition  of  26  mole  %  SiC. 

It  is  understandable  why  previous  attempts  to  codeposit 
HfB2  +  SiC  and  ZrB2  +  SiC  were  not  successful.  The  point 
identified  as  "A"  in  Figure  4  corresponds  to  the  point  which 
would  yield  HfB2  +  26  mole  %  SiC  if  the  deposition  efficiencies 
were  assumed  to  be  the  same  for  hafnium,  boron,  and  SiC.  Without 
benefit  of  the  thermodynamic  calculations,  this  reactant  composi¬ 
tion  would  be  a  logical  first  choice.  The  thermodynamically 
predicted  location  of  the  two  phase  HfB2  +  SiC  region  however,  is 
appreciably  shifted  to  the  right  (higher  MTS  concentration)  and 
down  (lower  BCI3  concentration)  compared  to  what  one  might 
anticipate.  Alternatively,  perhaps  the  deposition  temperature 
used  in  previous  experiments  was  too  high. 

We  have  also  performed  thermodynamic  calculations  for  the 
purpose  of  evaluating  the  influence  of  temperature  and  pressure 
on  the  codeposition  of  HfB2  +  SiC.  Briefly,  the  results  predict 
that  codeposition  of  the  desired  composite  should  be  possible 
over  the  temperature  range  1200-1500  K  at  pressures  from  0. 1-1.0 
atm.  Substitution  of  B2Hg  for  BCI3  appears  feasible.  Detailed 
results  of  these  calculations  were  presented  as  an  invited  paper 
at  the  Tenth  International  Conference  on  Chemical  Vapor  Deposi¬ 
tion.  [1]  More  recently,  additional  details  were  presonted. [ 13] 
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IV.  BASIC  KINETIC  AND  TRANSPORT  PRINCIPLES  AND  EQUATIONS 

The  purpose  of  modeling  chemical  vapor  deposition  is  to 
relate  operating  variables  (e.g.,  pressure,  temperature,  reactant 
concentrations,  reactor  geometry,  and  substrate  location  and 
orientation)  to  performance  parameters  (e.g.,  growth  rate, 
uniformity,  composition,  phase,  and  microstructure  of  the 
composites) . [14]  Such  models  can  be  constructed  by  describing 
the  underlying  principles  of  Kinetics  and  transport  phenomena  of 
CVD  processes  in  terms  of  solvable  mathematical  expressions. 
Generally  a  set  of  coupled-partial  (or  -ordinary)  differential 
equations  (i.e.,  momentum,  heat,  and  mass  balances)  is  required 
to  portray  a  CVD  process  realistically .[ 15 )  The  general  forms 
of  these  equations  which  are  applicable  to  most  CVD  processes  are 
introduced  in  this  section. 

Since  the  film  growth  rate  is  slew  compared  to  gas  phase 
dynamics  in  most  CVD  reactors,  steady  state  can  be  assumed.  Then 
the  equation  of  continuity  for  a  fluid  element  may  be  written: 

V-  (-^Tr)  =r  O  (1) 

where  V  *  vector  differential  operator 
a  total  density  of  the  fluid 
V  *  mass  average  velocity  of  the  fluid. 

The  momentum  balance  is: 

=  V*"  -  7  P  (2) 

where  -»  gravitational  acceleration 
p*  «  pressure. 

The  viscous  tensor  for  a  Newtonian  fluid  is: 

+  <.7ir)TJ  '  I3  I  <3> 

where  =  viscosity 

r  *  unit  matrix. 
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The  energy  balance  takes  the  form: 

“  V-(AVT)  +  HvRv) 


where  Cp  ■  specific  heat 
T  *  temperature 

3*  *  flux  of  species  i  relative  to  the  mass  average 
H*  *  partial  molal  enthalpy 

■  rate  of  production  of  species  i  due  to  homogeneous 
chemical  reactions. 

The  above  equation  states  that  convective  heat  transfer  is 
balanced  by  a  combination  of  heat  conduction  and  heat  generation 
by  inter-diffusion  of  species  and  gas  phase  reactions.  The 
Dufour  effect  and  viscous  energy  dissipation  are  neglected  since 
these  are  known  to  be  small  in  most  CVD  systems.  [16]  Heat 
evolved  from  heterogeneous  reactions  at  the  substrate  surface 
enters  into  Eq.  (4)  as  a  boundary  condition  (i.e.,  constant  heat 
flux  at  the  surface). 

The  steady  state  mass  balance  in  the  gas  phase  for  species 

i  is: 


where  N,  -  Jl  +  ^ 


(6) 


and  Cx.  is  the  concentration  of  species  i.  The  summation  of  these 
fluxes  is: 

£m,T, -o  (7) 

w*» 

where  M;  is  the  molecular  weight  of  species  i.  The  flux,  Nj,  can 
be  related  to  the  mole  fraction  of  species  i  inside  the  reactor 
by  (n-1)  Stefan-Maxwell  equations: 
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(8) 


r-  ft 

rrxi*Z_ - - - 

U  ^ 

where  Dv^  is  the  binary  diffusion  coefficient  between  species  i 
and  k  and: 


Ix\  =  ‘-0 


<*) 


The  boundary  conditions  for  Lqs.  (8)  and  (9)  are  the  surface 
fluxes  which  are: 


Mv 


Sur'ffcCft. 


r*»i 


(10) 


where  is  a  stoichiometric  coefficient  of  species  i  for 

heterogeneous  chemical  reactions ,  m,  with  the  reaction  rate  of 


Lm* 


The  rates  of  surface  reactions  have  the  general  form: 

y\ 


(11) 


where  is  a  Langmuir-Hingshelwood  expression.  This  has  the 

general  form  of: 


1«»>m 


(12) 


i  ^  Kwv.  C  Y'S) 

where  Kmt  is  the  equilibrium  constant  for  adsorption  and  desorp¬ 
tion  of  species  i  for  heterogeneous  reaction  m.  The  exponents  q 
and  n  are  determined  by  the  nature  of  the  reaction  mechanism. 

The  net  rate  of  production  of  each  species  in  gas  phase 
reactions  is  given  by: 


=  ~X  S-.«k 


(13) 


where  is  a  stoichiometric  coefficient  of  species  i  in  homogen¬ 
eous  reaction  n.  The  rate  of  homogeneous  reaction,  rn,  has  the 
general  form  of: 
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(14) 


where  V,  ■  forward  rate  constant  for  homogeneous  reaction  n 
■  backward  rate  constant  for  homogeneous  reaction  n 
Exact  forms  of  Eqs.  (10)  through  (14)  must  be  determined  (or 
assumed)  before  attempting  to  solve  Eqs.  (1)  through  (9).  The 
complexity  of  Eqs.  (10)  through  (14)  is  governed  by  kinetics. 

Eqs.  (1)  through  (4)  become  ordinary  differential  equations 
(ODE)  for  one-dimensional  analysis  and  partial  differential 
equations  (PDE)  for  multi-dimensional  problems.  It  is  emphasized 
that  solving  coupled-nonlinear-PDE  systems  is  a  difficult  task, 
if  at  all  possible. 

To  date,  two  "sophisticated"  models  have  been  developed 
using  this  general  approach  (i.e.,  simultaneous  solution  of  the 
above  equations)  after  considerable  mathematical  simplifications. 
Rebenne  and  Pollard[14]  have  used  a  finite  difference  method  to 
solve  momentum,  heat,  and  mass  equations  derived  for  CVD  of  boron 
carbide  in  an  impinging  jet  reactor.  Relatively  simple  kinetics 
are  assumed  for  the  above  modeling  purpose.  The  unique  features 
of  their  work  are  that  the  model:  (i)  considers  deposition  of 
mixed  solid  phases  (e.g.,  composition  and  relative  amount  of 
carbon  and  boron-carbon  phases)  and  (ii)  includes  the  largest 
number  of  chemical  species  (e.g.,  6  species;  BCI3,  CH4,  H2, 

BHCI2,  BC1 ,  and  HC1)  compared  to  other  CVD  models  to  date. 

Kinetic  parameters  such .as  rate  constants  are  estimated  by 
fitting  experimental  data  available  in  the  literature. 

Roenigk  and  Jensen  [17,18]  developed  a  2-dimensional  model 
for  low  pressure  chemical  vapor  deposition  (LPCVD)  of  silicon  and 
silicon  nitride  in  a  "hot-wall  multiple-waf ers-in-tube"  reactor. 
This  model  predicts  the  effects  of  operating  conditions  on  the 
deposition  rate  at  various  substrate  positions.  A  numerical 
method,  orthogonal  collocation,  is  used  to  solve  two  governing 
sets  of  ordinary  differential  equations  (e.g.,  one  for  the 
annulus  region  and  the  other  for  the  space  between  the  wafers). 
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V.  CHEMISTRY  OF  CVD  OF  BORON  NITRIDE  AND  ALUMINUM  NITRIDE 

As  mentioned  earlier,  in  order  to  model  a  multi-reaction 
system,  it  is  necessary  to  express  each  individual  reaction  step 
in  the  forms  of  Eqs.  (11)  and  (12)  for  heterogeneous  reactions 
and  Eqs.  (13)  and  (14)  for  homogeneous  reactions.  For  example, 
Roenigk  and  Jensen[17]  have  assumed  the  following  simple  kinetic 
expression  for  LPCVD  of  silicon  from  silane: 

- - - — 

This  expression  corresponds  to  a  reaction  mechanism  of:  (i) 
adsorption/desorption  equilibria  of  silane  and  hydrogen,  (ii)  no 
homogeneous  reactions,  and  (iii)  rate  controlling  surface 
reaction  (i.e.,  surface  decomposition  of  silane).  Kinetic 
parameters  k,  KH,  and  Kg  are  estimated  by  fitting  experimental 
data  using  nonlinear  regression.  The  above  expression  is  only 
valid  at  low  pressures  (ca.  10~2  Pa)  since  it  is  shown  that  gas 
phase  reactions  (e.g.,  pyrolysis  of  SiH4)  take  place  at  higher 
pressures  (19,20], 

CVD  of  silicon  is  a  much  simpler  and  better  understood 
process  than  that  of  boron  nitride  (BN)  or  aluminum  nitride 
(AIN) .  However,  even  for  this  relatively  simple  SiH4-Si-H2 
system,  a  mechanism  containing  20  intermediate  reaction  steps  has 
been  proposed  [21].  Therefore,  it  is  expected  that  one  of  the 
most  difficult  tasks  of  modeling  CVD  of  BN,  AIN,  and  BN  +  AIN  is 
understanding  the  chemistry  of  these  processes  and  developing 
(postulating)  possible  reaction  mechanisms.  It  seems  these 
mechanisms  need  to  be  deduced  based  on  our  own  experimental  data 
since  kinetics  of  BN-and  A1N-CVD  have  not  been  thoroughly 
reported  in  the  literature. 

However,  our  preliminary  experimental  results  have  already 
provided  new  insights  into  the  adsorption  behavior  of  BCI3  and 
AICI3  on  various  substrates.  It  is  observed  that  BN  and  AIN  can 
not  be  deposited  on  graphite  at  1400  K,  but  they  can  be  deposited 
on  alumina  and  molybdenum  at  1400  K.  This  behavior  is  explained 
in  Figure  6.  It  is  not  yet  sure  whether  homogeneous  reactions 
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between  BCI3  and  NH3,  which  produce  stable  intermediate  species 
such  as  CI2B-NH2  and  borazine,  take  place  in  the  gas  phase. 
Nevertheless,  it  is  suspected  that  a  filled  electron  orbital 
(i.e.,  electron  pair  donor)  existing  on  nitrogen  atoms  of  either 
NH3  or  the  above  mentioned  incermediate  species  favors  the  Lewis 
acidic  sites  (i.e.,  electron  pair  acceptor)  on  alumina  as 
preferential  adsorption  sites  as  shown  in  Figure  6a.  Once  a 
monolayer  of  BN  is  deposited  on  the  substrate  surface,  a  faster 
deposition  process  via  nitrogen-to-boron  and  boron-to-nitrogen 
adsorptions  is  expected  to  proceed  as  shown  in  Figure  6b.  This 
postulate  is  supported  by  the  fact  that  the  hexagonal  BN  layers 
are  packed  directly  on  top  of  each  other  and  a  boron  atom  in  one 
layer  of  these  hexagons  is  located  over  a  nitrogen  atom  in  the 
next  layer. [22]  This  catalytically  enhanced  behavior  is  consis¬ 
tent  with  experimental  observations  made  by  Takahashi  et  al.[23] 
They  have  shown  that  the  deposition  rate  of  BN  increases  by  a 
factor  of  10  when  carbon  steel  substrates  are  used  in  the  place 
of  alumina  due  to  the  presence  of  catalytically  active  metallic 
sites.  Obviously,  graphite  does  not  contain  these  acidic  sites 
(see  Figure  6c)  which  are  required  to  facilitate  surface  reac¬ 
tions  (rearrangements  or  recombinations)  into  the  BN  hexagonal 
structure.  * 

It  is  interesting  to  note  that  BN  can,  however,  be  deposited 
on  graphite  at  higher  temperatures  (ca.  2000  K).[24]  Therefore, 
it  seems  that  the  above  described  adsorption  mechanism  is  only 
valid  at  temperatures  around  1400  K.  At  2000  K,  thermal  energy 
may  be  large  enough  to  alter  the  above  mentioned  electronic 
interactions , 

VI.  MODELS  FOR  OUR  CVD  PROCESSES 

The  location  and  orientation  of  a  flat  substrate  to  be 
coated  in  our  CVD  reactor  are  schematically  shown  in  Figure  7. 

In  this  experimental  arrangement,  a  mixture  of  reactant  gases 
flows  into  the  furnace  cylinder  (reaction  chamber)  and  creates 
the  flow  pattern  (i.e.,  momentum  boundary  layer)  near  the  leading 
edge  of  the  flat  substrate  in  the  fluid  stream.  In  addition  to 
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the  momentum  boundary,  a  thermal  boundary  for  the  gas  mixture  and 
a  concentration  boundary  for  each  gas  species  develop  due  to  the 
presence  of  homogeneous  and  heterogeneous  reactions.  Then  the 
objective  of  modeling  is  to  determine  what  happens  when  heat, 
mass,  and  momentum  are  simultaneously  transferred  across  the 
boundary  of  the  flowing  fluid. 

Two  relatively  simple  (but  realistic)  models  appear  to  be 
promising  candidates  for  this  CVD  configuration.  They  are  the 
"boundary  layer  model"  and  the  "stagnant  layer  model".  Both  are 
based  on  the  boundary  layer  theory. [15]  The  stagnant  layer 
model  is  actually  a  simplified  version  of  the  boundary  layer 
model.  A  set  of  ordinary  differential  equations  is  required  for 
the  stagnant  layer  model  and  partial  differential  equations  are 
necessary  to  describe  the  boundary  layer  model.  Details  of  these 
models  are  discussed  in  this  section. 

1.  Boundary  Laver  Model 

For  the  boundary  layer  shown  in  Figure  7,  Eqs.  (1)  through 
(4)  are  simplified  to: 
continuity: 

m  +  = 

irx 

momentum : 

+- 

&  * 


energy: 


And,  Eqs.  (6)  through  (14)  remain  the  same.  The  boundary 
conditions,  if  not  stated  in  Eqs.  (6)  through  (13),  are: 


14' 


+  't-o  ■> 


~  O. 


(19) 


V a 

cd-  <^bo  x  £=.0  )  *Vi 


T-T~ 


x*. 


Xt 


uD 


(20) 

(21) 

(22) 

(23) 

(24) 


The  p^eceeding  equations  are  most  general  in  describing  our  CVD 
system.  More  simplifications  can  be  made  by:  (i)  assuming  the 
physical  properties  (e.g.,  ^  }  ve+o .  )  as  constants  if  the 

thermal  gradient  is  small  and  (ii)  treating  ”0^  as  small  compared 
to  'VJ  (i.e.  ,“VJ.  *0)  inside  of  the  boundary  layer.  For  a  hot  wall 
reactor  as  ours,  no  heat  flux  (i.e.,  no  energy  balance  and 
constant  physical  properties)  is  often  a  good  assumption  as  long 
as  heat  produced  or  absorbed  from  chemical  reactions  is  small. 

2.  Stagnant  Laver  Model 

Even  after  all  the  simplifications,  we  still  have  a  formid¬ 
able  task  of  solving  the  coupled  PDE  (e.g.,  Eqs .  (15),  (16), 

(18),  and  the  Stef an-Maxwell  equations).  This  difficulty  can  be 
avoided  by  introducing  the  stagnant  layer  concept.  This  assumes 
that:  (i)  reactants  are  not  depleted  along  the  axial  direction 

(i.e.,  constant  deposition  rate  along  the  axial  direction)  and 
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(ii)  velocity,  thermal,  and  concentration  profiles  are  the  same 
at  any  axial  position  (i.e.,  no  variation  in  tho  boundary  layers 
thickness),  as  shown  in  Figure  8.  This  model,  therefore,  can  not 
predict  the  deposition  rate  as  a  function  of  the  axial  direction 
which  may  be  an  important  design  consideration  in  industrial 
applications.  However,  it  is  believed  that  the  length  (ca.  one 
to  three  inches)  of  the  substrates  used  in  our  experiments  is  too 
short  for  detecting  large  variations  in  film  thickness. 

The  governing  equations  then  are: 


energy: 


mass : 


SR, 

d* 


7  ( T  *iik 


RR»)  = 


=  o 


(25) 


(26) 


with  Eqs,  (6)  through  (14). 

Using  this  stagnant  layer  model  concept,  a  very  idealised 
situation  for  CVD  of  BN  from  BCI3-NH3-H2  mixtures  at  typical 
operating  conditions  was  modeled.  The  following  assumptions  were 
made: 

(i)  no  thermal  energy  gradient, 

(ii)  mass  transfer  through  the  stagnant  layer  by  molecular 
diffusion  only, 

(iii)  constant  physical  properties, 

(iv)  decomposition  of  BCI3  at  the  surface  is  the  rate  control¬ 
ling  step  so  the  mass  transport  and  surface  decomposition 
of  BCI3  are  two  competing  processes  which  determine  the, 
overall  deposition  rate,  and 

(v)  no  gas  phase  reactions. 

With  the  above  assumptions,  substitution  of  Eq.  (6)  into  Eq.  (26) 
gives : 

.  Ok 


X  = 


-'ae  f 
RT  V 


-) 


(27) 
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where  is  the  diffusion  coefficient  of  BCI3  in  the  fluid 

mixture.  The  first  order  approximation  of  the  above  equation  is: 

X*  <28> 

where  f^u  and  are  the  partial  pressures  of  BC13  in  the 

bulk  phase  and  at  the  surface,  respectively,  and  is  the 
thickness  of  the  stagnant  layer  of  BCI3.  From  assumption  (iv), 
the  rate  controlling  step  is: 

t.  •  ri. 

The  above  kinetic  expression  is  derived  based  on  the  pseudo  first 
order  surface  decomposition.  At  a  steady  srate,  the  overall 
growth  rate, 3*^  •  is  equal  to  andT^: 

X  -  x  -  X  <3°> 


Combining  Eqs.  (28),  (29)  and  (30)  gives: 

XT  =  £§I -  (31) 

VDab  +  '/^ 

The  terms  V /5>ab  and  */"&.  can  de  thou9ht:  as  "diffusion  resis¬ 

tance"  and  "kinetic  resistance"  to  the  process,  respectively.  In 
other  words,  the  process  is  limited  by  diffusion  if  muc^ 

larger  than  (/-(^  and  is  kinetically  controlled  if  the  opposite  is 
true.  This  simple  concept  was  first  introduced  by  van  den 
Brekel.[25]  He  combined  the  above  two  terms  and  defined  a 
dimensionless  number  called  the  CVD  number: 


CVD 


iv 

Dab 


(32) 


Using  the  above  equations  and  the  physical  properties  obtained 
from  the  thermodynamic  correlations  given  in  Reid  et  al.[26], 
the  CVD  number  for  CVD  of  BN  is  calculated  to  be  about  8.5  at  the 
following  conditions: 
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P  -  1  atm,  T  -  1400  K,  total  flow  rate  -  650  SCCM 
PNH3  -  0.65  atm,  PBC13  «  0.04  atn,  and  PH2  -  0.31  atm 
The  rate  constant  k  is  estimated  from  experimental  work  of  Male 
and  Salanoubat. [27]  This  result  (i.e.,  CVD  number  ■  8.5) 
suggests  that  the  overall  deposition  rate  is  limited  by  both 
diffusion  and  kinetic  steps  while  the  effect  of  diffusion- 
resistance  is  greater.  As  soon  as  more  experimental  data  are 
available,  calculations  will  be  extended  to  account  for  multi- 
component  diffusion  and  multi-homogeneous  and  -heterogeneous 
reactions. 

3.  Heat  Transfer 

As  mentioned  earlier,  the  presence  of  the  thermal  boundary 
(or  stagnant  layer)  may  play  an  important  role  in  determining  the 
complexity  of  mathematics .  A  preliminary  calculation  was  per¬ 
formed  to  discern  whether  or  not  the  “heat  transfer  effect"  is 
important  in  CVD  of  BN  and  AIN.  For  this  purpose,  we  have 
assumed  that  reactions  take  place  at  the  surface  only.  The  heats 
of  reaction  data  were  obtained  from  the  Janaf  Thermochemical 
Tables. [28] 

BC13  +  NH3  ->  BN(s)  +  3HC1) AHrxn  *  -17.8  kcal/mol  at  1400  K 
A1C13  +  NH3  -*  AlN(s)  +  3HCi)4Hrxn  *  10.5  kcal/mol  at  1400  K 

As  can  be  seen  from  the  above  equations,  CVD  of  BN  is  an  exother¬ 
mic  process  so  heat  is  generated  at  the  surface  as  the  reaction 
proceeds  (i.e.,  constant  heat  flux  at  the  surface).  The  equation 
governing  the  laminar  heat  transfer  along  a  flat  plate  which  is 
at  constant  heat  flux  is  given  by[29]: 

*  average  temperature  difference  between  the 
plate  and  bulk  phase  over  the  plate  length,  L. 
ReL  =  Reynold's  number  at  the  length  L 
Pr  =  Prantle's  number 

qw  *  heat  flux  (&Hrxn  x  deposition  rate) 


where 


From  the  Above  equation,  the  average  temperature  difference 
is  estimated  to  be  about  32  K  at  a  typical  operating  condition 
(ca.  growth  rate  ■  1  m/min,  1400  K,  1  atm,  total  flow  rate  -  450 
SCCM,  XBC1  *  0.056  atm,  ■  0.222  atm,  «  0.056  atm,  and 

xAr  ■  0.66^  atm).  This  suggests  that  the  effect  of  heat  transfer 
may  not  be  neglected  in  future  studies.  However,  the  effect  of 
heat  transfer  for  CVD  of  AIK  is  expected  to  be  much  smaller  than 
that  of  BN  because  of  the  low  AIN  deposition  rate. 

VII.  FUTURE  MODELING  EFFORTS 

As  mentioned  previously,  reasonable  kinetic  mechanisms  are 
required  for  developing  a  complete  mathematical  model.  However, 
the  complexity  of  CVD  processes  for  BN,  AIN,  and  BN  +  AlN  is  the 
biggest  obstacle  in  determining  tho  reaction  mechanism  associated 
with  these  processes.  Therefore,  a  major  portion  of  future 
experiments  will  be  devoted  to  understanding  th*  chemistry  of 
these  processes  (i.e.,  mainly  the  effects  of  reactant  concentra¬ 
tions,  temperature,  and  pressure  on  performance  parameters). 

Also,  a  numerical  algorithm  for  solving  a  coupled-ODE  system  will 
be  constructed  and  finalized  while  the  above  discussed  experi¬ 
ments  are  concurrently  conducted. 

VI I I.  COATING  EQUIPMENT 

An  automated  system  for  introducing  reagents  into  the 
Georgia  Tech  CVD  System  was  procured.  This  equipment,  which 
consists  of  four  mass  flow  meters  and  a  pressure  control  system, 
permits  accurate  control  of  reagent  flow  rates;  and  gas  pressure. 
With  these  modifications  the  coating  system  is  fully  automated 
and  capability  exists  for  conducting  controlled  experimentation. 

A  mass  spectrometer  will  soon  be  purchased  which  will  permit 
campling  and  analysis  of  gases  within  the  coatirg  reaction 
chamber.  This  will  provide  qualitative  and  quantitative  data  on 
gaseous  species  present  within  the  reactor.  This  information 
will  greatly  assist  CVD  model  development  as  well  as  permit 
independent  monitoring  of  the  performance  of  the  reagent  supply 
system. 


19 


IX.  PRELIMINARY  EXPERIMENTAL  RESULTS 

Fourteen  CVD  experiments  were  performed  with  the  intention 
of  depositing  BN  or  AIN  or  codepositing  BN  +  AIN.  The  substrate 
materials  employed  in  these  experiments  included  various  shapes 
of  graphite,  alumina,  fused  silica,  and  silicon  nitride.  The 
substrates  were  suspended  by  molybdenum  wire.  The  temperature  in 
the  reaction  chamber  was  1400  R  for  each  run.  Run  condition 
variables  included  the  rates  of  reagent  gas  flow,  total  gas 
pressure  within  the  reaction  chamber,  and  the  duration  of  the 
run.  Some  deposition  resulted  in  the  majority  of  these  experi¬ 
ments.  Examination  by  light  microscopy  and  ESCA  is  complete  but 
additional  characterization  is  in  progress.  Most  of  the  samples 
will  be  studied  by  scanning  electron  microscopy,  and  selected 
samples  will  be  analyzed  by  XRD  and  the  electron  microprobe.  The 
following  paragraphs  are  a  summary  of  visual  observations  and 
ESCA  results. 

Four  BN  deposition  experiments  were  conducted  using  BCI3, 
NH3,  H2,  and  Ar  as  reagents.  At  one  atmosphere  pressure  a 
white,  uniform  coating  was  obtained  on  alumina  and  silicon 
nitride  substrates  and  on  the  molybdenum  wire.  An  increase  in 
the  NH3/BCI3  ratio  caused  little,  if  any,  change  in  the  appear¬ 
ance  of  the  deposit.  The  deposit  on  alumina  and  silicon  nitride 
became  thicker  and  more  translucent  in  an  identical  period  of  run 
time  when  the  pressure  in  the  reaction  chamber  was  reduced  to  0.1 
atm.  One  attempt  to  deposit  BN  using  the  reagents  BCl3,H2/  and 
N2  was  made  but  no  coating  was  obtained.  This  indicates  that  N2 
cannot  be  substituted  for  NH3  as  a  nitrogen  source. 

The  first  two  of  five  AIN  deposition  experiments  produced  no 
coating.  The  first  run  was  aborted  due  to  clogging  of  the  gas 
injector.  The  other  experiment  which  produced  no  coating  may  not 
have  had  a  sufficient  amount  of  A1C13  vaporized.  The  temperature 
of  the  vaporizer  for  that  experiment  was  146°C.  The  three 
subsequent  experiments  resulted  in  deposits.  The  temperature  of 
the  AICI3  vaporizer  was  increased  to  152°C  with  a  run  time  of  60 
minutes  for  two  of  these  deposition  runs.  The  best  coating 
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resulted  when  the  temperature  of  the  AICI3  was  increased  to  162*C 
with  a  run  time  of  40  minutes 

The  first  of  four  BN  +  AIN  codeposition  experiments  was 
aborted  after  20  minutes  due  to  clogging  of  the  gas  injector.  A 
small  amount  of  white  material  was  deposited  on  the  substrates. 
The  ammonia  flow  rate  was  reduced  and  the  argon  flow  was  greatly 
increased  in  the  next  experiment,  which  resulted  in  a  non- 
uniform,  beige  colored  deposit.  "Mushroom"  shaped  protrusions  on 
the  coating  were  observed  using  an  optical  microscope.  A  more 
uniform,  beige  colored  deposit  which  was  crystalline  in  appear¬ 
ance  resulted  when  hydrogen  was  added  to  the  reagent  gas  stream 
in  the  next  experiment.  For  the  last  experiment,  the  flow  rates 
of  ammonia  and  hydrogen  were  increased  and  thin  beige  colored 
deposits  resulted.  Examination  of  these  coatings  by  ESCA 
revealed  about  the  proper  ratios  for  the  B,  Al,  and  N  contents. 
Additional  characterization  will  emphasize  determination  of  the 
phases  present  and  coating  microstructure. 
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Deposition  diagram  for  the  reactant  system  BCI3-AICI3-NH3  showing 
e  BN  +  AIN  can  be  codeposited. 


■v*  w  rw/w/»wwnfw^> 


ivdnvrvirj  Tv;ritiru:iu.r 


Figure  3.  Partial  pressures  of  the  major  gaseous  species  as  a 
function  of  reactor  temperature  for  the  reactant  system  BCI3- 
AICI3-NH3  which  is  used  for  codeposition  of  BN  +  AIN.  The 
relatively  high  stability  of  the  A1C1X  species  compared  to 
species  containing  boron  accounts  for  the  deposition  efficiency 
of  AIN  being  less  than  that  for  BN. 
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Figure  4 .  CVD  deposition  diagram  for  the  reactant  system  HfCl4 
BCI3-MTS  at  1200  K  showing  region  where  H£B2  +  Sic  codeposit. 
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Figure  5,  CVD  deposition  diagram  for  the  reactant  system  HfCl4 
BCI3-MTS  at  1200  K  showing  mole  %  SiC  contour  lines  in  the  two 
phase  region  HfB2  +  SiC. 
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Figure  6.  Adsorption  behavior  of  boron  containing  gas  species  on  various  sub¬ 
strates  (a)  alumina  substrate,  (b)  BN  substrate,  (c)  graphite  substrate. 
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Figure  7.  The  laminar  momentum  boundary  layer  along  a  flat 
substrate  located  inside  the  CVD  reactor. 


